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Abstract. Recent studies have reported declines in lizard populations associated with local warming trends and alteration 
of thermal niches. These studies suggest that there are some key physiological and behavioural parameters that determine 
the sensitivity of each species to the local context of global warming. The Atlas Dwarf Lizard, Atlantolacerta andreanskyi, 
is a small lacertid lizard endemic to the High Atlas Mountains (2200-3500 m a.s.l.) of Morocco. Its populations display a 
disjointed distribution across 440 km, consisting of phylogenetic lineages separated since the Miocene, suggesting a mor- 
phologically cryptic species complex. To establish whether this deeply rooted evolutionary divergence has promoted shifts 
in fundamental niche occupancy, but also whether there is any ecophysiological plasticity, preferred body temperatures 
(Tpref) and evaporative water loss (EWL) were examined in populations of the central and western High Atlas. Our results 
revealed no differences in Tpref across populations, but lower thermal preferences in (gravid) females. In contrast, EWL 
varied between populations, but not between sexes. In both cases, sexual patterns remained similar between populations, 
and differences were not evident between body sizes. Importantly, Tpref and EWL traded off, even if with a common trend 
across populations. These results suggest that ecophysiology remains mostly conservative within this species complex and 
that hydric ecology, which is more closely associated with body size, is more flexible than thermal ecology. The increase of 
lizard body size with altitude, in agreement with life history trends with seasonality observed in other Mountain Lizards, 
cannot be interpreted in terms of ecophysiological adaptation since the low-altitude populations of small lizards are the 
ones most exposed to dehydration. Rather, local environment and possibly aspects of evolutionary history could be in- 
volved. Overall, the ecophysiological conservativeness of this endemic species confirms its vulnerability to climate change, 
and indicates that hydric ecology needs to be incorporated into general models regarding how behavioural plasticity may 
mediate resilience to such changes. In this context, the low-altitude populations are the most vulnerable ones, and local 
management measures should focus on the conservation of vegetation and water bodies. 
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Introduction 


Global climate change and associated alterations to ther- 
mal niches are predicted to have a catastrophic impact on 
lizard diversity, with models suggesting local extinctions 
reaching 39% being validated by observed local extinction 
rates (SINERVO et al. 2010). It is therefore essential to estab- 
lish how plasticity of thermal ecology within lizard species 
may mediate their resilience to climate change. Temper- 
ature and humidity are among the most important envi- 
ronmental factors that define a reptile’s niche, influencing 


metabolic reactions, heat and water balance within a giv- 
en environment, and ultimately the performance and fit- 
ness of individuals (PORTER & TRACY 1983, ANGILLETTA 
2009, SANNOLO & CARRETERO 2019). Unlike endotherms, 
which produce heat to maintain their body temperature, 
ectotherms, including lizards, typically thermoregulate 
through behavioural adjustments (STEVENSON 1985, Coss- 
INS & BOWLER 1987). For lizards, it was hypothesized that 
the preferred body temperature (Tpref) would be either an 
evolutionary labile or static trait (HERTZ et al. 1983). The 
labile hypothesis supposes fast local adaptation to the ther- 


© 2022 Deutsche Gesellschaft für Herpetologie und Terrarienkunde e.V. (DGHT), Germany 


Open access at https://www.salamandra-journal.com 


139 


ABDERRAHIM S’KHIFA et al. 


mal environment, while the static hypothesis emphasises 
the evolutionary rigidity of the Tpref, with adjustments to 
local conditions being mediated by phenotypic plasticity, 
and in particular by microhabitat selection or activity pat- 
terns (GvozDik 2002). For a long time, most studies on liz- 
ards were focused exclusively on the thermal aspect of their 
ecology, but recent studies have since illustrated its interde- 
pendency with water balance; for example demonstrating 
that dehydrated lizards tend to prefer lower body tempera- 
tures (CROWLEY 1987, SANNOLO & CARRETERO 2019). 

In lacertid lizards, Tpref varies in response to season, 
reproductive condition, and feeding status (CASTILLA et al. 
1999, CARRETERO et al. 2005). However, within compar- 
ative classes, lacertid Tpref is generally considered evolu- 
tionarily rigid (Van Damm et al. 1990, but see SAGONAS 
et al. 2013), responding very slowly to long-term changes 
in the thermal environment (GARCIA-PORTA et al. 2019). 
As a consequence, Tpref tends to remain similar within 
species even under different climate regimes (VAN DAMME 
et al. 1989, 1990, TOSINI et al. 1993, GVOŽDÍK & CASTILLA 
2001, CARRETERO et al. 2005). This is because behavioural 
thermoregulation is expected to buffer the effects of ther- 
mal heterogeneity in the environmental context, restrain- 
ing adaptation in thermal physiology (LoGAN et al. 2019). 
Nevertheless, recent studies found individual variation in 
thermal preferences within species, suggesting that natu- 
ral selection may still act on these traits (ARTACHO et al. 
2013, 2015, MUÑOZ & BODENSTEINER 2019). Consequently, 
both interspecific and intraspecific variation in Tpref must 
be analysed taking into account the phylogeny and biogeo- 
graphic history of the groups under study. Different ther- 
mal responses depend on variable thermal requirements 
due to ontogenic state, sex or reproductive cycle (STEVEN- 
SON 1985, ATKINSON 1994, CULLUM 1998, CLUSELLA-TRUL- 
LAS et al. 2009). In particular, reproductive females often 
differ in their thermoregulatory behaviour, since embryon- 
ic development requires a particular thermal regime that 
will often differ from that of adults (GEORGES et al. 2005). 

Adding an extra layer of complexity, many lizard spe- 
cies range across steep altitudinal gradients. Generally, 
environmental temperature decreases as altitude and lat- 
itude increase, which makes the thermal environment at 
high latitudes and altitudes (e.g., high-mountain systems) 
particularly hostile to ectotherms (Gvozpixk 2002, ANGIL- 
LETTA 2009, AGUADO & BRANA 2014). In such cold envi- 
ronments, ectotherms may not be able to thermoregulate 
effectively for long periods and, hence, they must either 
be efficient thermoregulators or become inactive (HUEY & 
SLATKIN 1976, VAN DAMME et al. 1989, HERTZ et al. 1993, 
ANGILLETTA et al. 2002, ANGILLETTA 2009). Since ther- 
mal physiology may respond slowly to such pressures and 
thermal requirements remain conservative in populations 
and species (HERTZ et al. 1983, CROWLEY 1985, ANDREWS 
1998, SEEBACHER & SHINE 2004, DIAZ DE LA VEGA-PEREZ 
et al. 2013), species inhabiting contrasting environments 
(e.g., high vs. low altitudes) often exhibit similar thermal 
requirements but incur different thermoregulatory costs 
(BOGERT 1949, HERTZ et al. 1983, CROWLEY 1985, S KHIFA 
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et al. 2020a). In such cases, the organisms must compen- 
sate for poor thermal quality by highly efficient thermoreg- 
ulatory behaviour or other phenotypic shifts (DIAZ 1997, 
ASHTON & FELDMAN 2003, OSOJNIK et al. 2013). 

Thermoregulating species can respond to thermal het- 
erogeneity by selecting specific microhabitats and body 
postures to optimise heat gain and loss (HUEY & BENNETT 
1987). While doing so, they will also experience water loss 
through evaporation, and since evaporation increases with 
rising mean temperatures, a compromise is made between 
thermoregulation and water balance (MauTz 1982a, b). 
Thus, dehydration could pose another constraint to activity 
when the availability of water is limited. Indeed, recent re- 
search on lizards suggests that humidity, ambient moisture, 
and the availability of water might be even more limiting 
than the thermal environment in warm areas (CARRETERO 
et al. 2016, FERREIRA et al. 2016, S KHIFA et al. 2020a). Also, 
while their thermal preferences may be similar, lizard pop- 
ulations with limited access to water may be more resist- 
ant to water loss than those populations with unrestricted 
access to water (DUPOUE et al. 2017). Moreover, gravidity 
may also affect water loss because embryo needs may shift 
the balance between thermoregulation and hydroregula- 
tion in the mothers. While ecophysiological shifts in gravid 
females were traditionally interpreted as being due to the 
higher oxygen needs of large embryos at higher tempera- 
tures (MADERSON & BELLAIRS 1962), there is now evidence 
that dilated bellies may increase water loss rates, inducing 
mothers to counteract the dehydration of their offspring 
(Bonner et al. 2017), with both probably acting in combi- 
nation. However, generalisation of these findings is prob- 
lematic due to a lack of comparative data, since most of the 
current literature on lizard ecophysiology remains focused 
on thermal traits (HERRANDO-PEREZz et al. 2020). For ex- 
ample, it is almost unknown to which extent water loss 
rates may vary along environmental gradients, or whether 
they are subject to diel and seasonal plasticity, as thermal 
preferences often are (SANNOLO et al. 2020). Incorpora- 
tion of such information into modelling approaches will 
be essential to establish whether behavioural responses of 
lizards may mediate the worst impacts of global climate 
change. 

Here, we investigated the thermal and hydric ecophysio- 
logy of the Atlas Dwarf Lizard Atlantolacerta andreanskyi 
in contrasting mountain systems of the High Atlas range. 
We compared the Tpref and evaporative water loss (EWL) 
rates between three populations representing widely di- 
vergent lineages within this species complex (BARATA et 
al. 2012). In accordance with recent findings in the fam- 
ily Lacertidae (GARCIA-PORTA et al. 2019, SANNOLO et al. 
2020), we would expect more flexibility across lineages in 
hydric physiology than in thermal physiology. Moreover, 
we predict that gravidity will induce phenotypic plasticity 
in females, affecting the balance between thermoregulation 
and hydroregulation. Furthermore, as the combined study 
of thermal and hydric ecology is the basis of the delimita- 
tion of fundamental niche occupancy in lizards (BOWKER 
1993), our aim was to establish whether temperature and 
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water requirements trade off, and if this is the case, if the 
same balance is maintained across populations/lineages. 
This information may be critical for validating correla- 
tive ecological models based on realized niche occupancy, 
which already indicate that precipitation plays a promi- 
nent role in the distribution of Moroccan endemic reptiles 
and forecast considerable range retractions due to climate 
change (MARTÍNEZ-FEIRÍA et al., 2013). 


Material and methods 
Species and areas studied 


The Atlas Dwarf Lizard, Atlantolacerta andreanskyi (WER- 
NER 1929) is a small lacertid lizard endemic to Morocco, 
where it is found from 2200 to 3500 m a.s.l. in the moun- 
tains of the High Atlas and the Siroua Massif (MARTINEZ 
DEL MARMOL et al. 2019). Its distribution (Fig. 1) is clearly 
discontinuous, and extends from Jbel Ayachi in the east to 
Jbel Aoulime in the west for about 440 km (straight line). 
Associated with this spatial disjunction, almost no cur- 
rent gene flow between the different populations has been 
identified (BARATA et al. 2012). A mountain specialist, this 
ground-dwelling lizard is often found near watercourses 
and around the bases of cushion-like thorny plants (Bons 
& GENIEZ 1996, SLIMANI et al. 1996) that offer a buffered 
microclimate with higher humidity, more food, and pro- 
tection against predators and wind (SCHLEICH et al. 1996). 
Genetic studies conducted by BARATA et al. (2012) and 
GIOVANNOTTI et al. (2020) suggest that A. andreanskyi 
may be a species complex with at least seven distinct line- 
ages grouped in two major clades, i.e., a western (Oukai- 
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meden + Jbel Siroua) and an eastern clade (Jbel Ayachi + 
Outabati + Jbel Azourki + Tizi-N-Tichka). Although the 
cladistic relationships of the Jbel Aoulime population are 
poorly supported, it may be the sister taxon to the rest of 
the eastern lineage (BARATA et al. 2012). These distinct lin- 
eages separated during the Miocene (7.6 Mya), while the 
separation between the seven lineages (likely species) in- 
habiting separated mountain massifs is estimated to be 
Plio-Pleistocenic. Contrasting with this deep phylogenetic 
divergence, morphological variation showed considerable 
overlap across lineages, limited diagnostic traits, and only 
some size differences, suggesting niche conservatism with- 
in this species complex (BARATA et al. 2015). 

In an ecophysiological study conducted in Oukaime- 
den, ŞKHIFA et al. (2020a) reported that Aflantolacerta 
had similar thermal requirements to sympatric species of 
the genera Podarcis and Scelarcis, but was more sensitive 
to dehydration. Given the long evolutionary time frame, 
it seems that aridification of climate, together with niche 
conservativeness due to the interaction with other lowland 
lacertids, may have resulted in its range retraction and 
fragmentation while preventing niche shifts. However, this 
hypothesis is based solely on evidence from Oukaimeden. 
Therefore, we here investigate if the fundamental niche oc- 
cupancy of A. andreanskyi has remained conserved in the 
long term by comparing representatives of three lineages. 
Some aspects of the biology of A. andreanskyi are already 
well known, at least at Oukaimeden (BUSACK 1987, CAR- 
RETERO et al. 2006), but it remains unclear how the differ- 
ent climatic events that occurred during the Miocene and 
Pleistocene may have shaped the ecophysiological varia- 
tion in this species complex. 
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Figure 1. Distribution map of Atlantolacerta andreanskyi (red dots) and localities studied (stars) in the High Atlas Mountains. 
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Populations representing these lineages (eastern, west- 
ern and the ambiguously placed Jbel Aoulime) were inves- 
tigated in three high-altitude localities (Fig. 1) in the High 
Atlas in Morocco. All are characterized by a mountain cli- 
mate with a cold winter and sport vegetation consisting 
of xerophilous thorny bushes (Alyssum spinosum) that of- 
fer excellent protection against predators: 1) Ait Bougue- 
mez (51 km southeast of Azilal in the Central High Atlas; 
31°4545.8” N, 6°16'50.6” W, altitude 2,723 m a.s.l. — this cor- 
responds to the population called “Jbel Azourki” in the pre- 
viously published genetic analyses), 2) Oukaimeden (75 km 
south of Marrakech on the southern limits of the Central 
High Atlas; 31°12’09.5” N, 7°5124.4” W, altitude 2,600 m 
a.s.l.), and 3) Jbel Aoulime (64 km south of Imintanoute 
in the Occidental High Atlas; 30°58’50.7” N, 8°45711.1” W, 
altitude 2,512 m a.s.l.). Climatic data (Supplementary Ta- 
ble S1) of the three study areas, gathered from www.power. 
larc.nasa.gov for the last 30 years (1990-2020), show that 
Jbel Aoulime was less humid and warmer than Oukaime- 
den and Ait Bouguemez during this period. The ombro- 
thermic diagrams (Supplementary Fig. S1) produced for 
the 30 years of recording for the three localities indicate a 
dry season (precipitation < 2*temperature) extending over 
5, 6, and 7 months, respectively, at Ait Bouguemez, Oukai- 
meden, and Jbel Aoulime. 


Sampling 


Sixty adult individuals (11 males and 9 females from Ait 
Bouguemez, 10 males and 10 females from Oukaime- 
den, and 8 males and 12 females from Jbel Aoulime) were 
captured by hand along random one-way transects dur- 
ing their activivity periods from 9:00 to 18:00 h GMT in 
the first week of June 2019. We recorded the snout—vent 
length (SVL) for each individual, using a digital calliper 
(precision 0.01 mm), and the body mass (BM) with a dig- 
ital scale (precision 0.01 g). Due to the sampling period, 
all females captured and used in our different experiments 
were gravid (apparent through palpation). After capture, 
the lizards were kept captive in individual plastic terraria 
(30 x 19 x 17 cm) for no more than five days with water and 
food (Tenebrio molitor larvae) available to them ad libitum 
daily during our ecophysiological experiments. The ex- 
periments were carried out in the Water, Biodiversity and 
Climate Change Laboratory, Biodiversity and Ecology of 
Vertebrates Unit, of the Faculty of Sciences Semlalia Mar- 
rakech, Morocco. 


Preferred body temperature 
and water loss experiments 


Evaporative water loss (EWL) and preferred body temper- 
ature (Tpref) were evaluated for a target sample size of 20 
adult lizards per population. Each lizard was subjected to 
the two experiments (once each). The experiments of water 
loss were always conducted after the experiments of pre- 
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ferred body temperature, with an interval of two days in 
between. For both experiments, the lizards were left to fast 
for a day before handling. All lizards were released at the 
exact locations of their capture after the end of the experi- 
ments. 

The experiments of Tpref were executed in a closed 
room in the laboratory to prevent wind and sun from af- 
fecting the temperature in the terrarium. Each lizard was 
placed in a thermal gradient corridor (min-max 13-50°C; 
100 X 30 X 30 cm) generated by a thermoSpotlight 100 W 
(HOBBY®), fixed 20 cm above the substrate at one extrem- 
ity of the gradient. The bulb was switched on one hour 
before a lizard was placed into the corridor, and the first 
measurement was taken one hour later. The experiments 
were conducted between 9:00 and 19:00 h GMT. Due to 
logistic constraints (no thermal camera available) and to 
prevent the disturbance of gravid females due to the per- 
manent or repeated introduction of thermal probes to the 
cloaca (Koziz1 et al. 2021), body temperature was meas- 
ured every hour on the dorsal skin of the lizards using an 
infrared thermometer (FLUKE 572 CF, precision 0.01°C), 
perpendicularly orientated to the lizard and at very close 
distances (approximately 10-20 cm) to the animal to en- 
sure accurate measurements of body temperatures (HARE 
et al. 2007) without handling the lizards. 

For the experiment of EWL, each lizard was housed in- 
dividually in a plastic box (15 x 10 x 3 cm, ~ 0.45 1) with 
ventilation holes on the top, bottom and sides. Each box 
was placed into a second box including 5 g of silica gel. All 
individual plastic boxes were placed in a large container 
(40 x 30 X 25 cm, ~ 30 1), with 100 g of silica gel in a gauze 
bag. The quantity of silica gel used facilitated a low level 
of relative humidity (20 ~ 30%) inside each box, where we 
installed a thermo-hygrometer (Beurer HM 16, precision 
0.1%; 0.1°C) to monitor humidity and temperature (Oso- 
JNIK et al. 2013, CARNEIRO et al. 2015, SKHIFA et al. 20204). 
The container with the boxes in it was then placed in a cli- 
matic chamber (POL-EKO-APARATUR SP. J type ST3+/ 
ST3+) set at 24°C. This temperature was selected to allow 
lizard activity while preventing acute desiccation stress, 
which may occur if lizards remain at Tpref for long peri- 
ods without an opportunity for thermoregulation (San- 
NOLO et al. 2018, SANNOLO & CARRETERO 2019). This ex- 
periment ran from 8:00 to 19:00 h GMT. Lizards were in- 
dividually removed from the container at hourly intervals, 
weighed together with the box using a digital scale (RAD- 
WAG, AS 110/C/2, precision + 0.0001 g), and immediately 
placed back into their respective containers in the climatic 
chamber. To minimize stress and hyperventilation, which 
may increase EWL, lizards’ weights were obtained by sub- 
tracting the weights of the respective individual boxes. The 
whole measuring procedure took no longer than 20 sec- 
onds for each lizard. 

From our data, we calculated three measures of relative 
water loss. The accumulated water loss (EWLa) after each 
hour of the experiment (from 8:00 to 19:00 h GMT), was 
calculated according to the formula: 

EWLa = [(Wo - Wn) / Wo] x 100, 
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Table 1. Sample size, mean, and standard deviation (SD) of the body mass (BM), snout-vent length (SVL), preferred body temperature 
(Tpref), total evaporative water loss (EWLt), and instantaneous evaporative water loss (EWLi) calculated for each sex and population 
of A. andreanskyi. Jbel Aoulime (JAO), Oukaimeden (OUK), and Ait Bouguemez (BOU). * Groups whose elements are homogeneous 


after performing post hoc tests at a = 0.05. 


oe BM (g) SVL (mm) Tpref (°C) EWLt (%) EWLi (%) 
E Localities Nümber mean + SD mean + SD mean + SD mean + SD mean + SD 
, Jbel Aoulime 12 1.83+0.69 46.50+2.61 33.2+1.1 3.41+1.02 0.31+0.09 
Sen Oukaimeden 10 1.86+0.27 47.90+2.23 34,5+0.9 4.23+41.77 0.38+0.16 
Ait Bouguemez 9 3.36+0.59 50.67+4.24 34.3+1.5 3.00+0.89 0.27+0.08 
Jbel Aoulime 8 1.66+0.49 43.38+4.43 34,5+1.7 4.16+1.96 0.38+0.18 
Males Oukaimeden 10 1.89+0.24 46.30+1.77 34,9+1.0 3.28+1.62 0.30+0.15 
Ait Bouguemez 11 2.65+0.27 49.18+2.48 35.0+1.0 2.29+0.65 0.21+0.06 
os P<0.001 P<0.001 P=0.016 P=0.016 
*Statistical group: P=0.055 
A {JAO-OUK} {JAO-OUK} {JAO-OUK} {JAO-OUK} 
between populations {BOU} {BOU} {JAO-OUK-BOU} {BOU} {BOU} 


where W is weight. The instantaneous water loss (EWLi) 
was calculated using the formula 

EWLi = [(Wn - Wn+1) / Wo] x 100. 

And the total water loss (EWLt) was calculated by sub- 
tracting the mass of the lizard at the end of the experiment, 
and dividing it by the initial mass at 8:00 h (Wo) and mul- 
tiplying it by 100 to obtain a percentage value (%). While 
EWLa provides information on the overall volume of wa- 
ter evaporated since the start of the experiment, EWLi val- 
ues record the pattern of water loss over time (constant or 
variable), and EWLt is the total amount of water lost by the 
animal at the end of the experiment (Osojnik et al. 2013, 
SKHIFA et al. 20204). 


Statistical analyses 


We first tested for morphological variation in body mass 
(BM) and snout-vent length (SVL) across populations and 
between sexes with a two-way ANOVA, as well as variation 
in robustness with the same factors by comparing body 
masses with SVL as a covariate in an ANCOVA. Morpho- 
logical variables were log-transformed for linearizing the 
relations. The residuals of the ecophysiological data were 
normally distributed (Shapiro-Wilk test, p > 0.05), homo- 
scedastic (multivariate Box M and univariate Levene’ tests, 
p > 0.05), and non-spherical (Mauchley’s tests, p > 0.05 in 
most cases). Two-way AN(C)OVA tests (populations and 
sex as between-subject factors) of Tpref and water loss total 
rate were performed based on the values of both ecophysio- 
logical parameters of males and females of each popula- 
tion. Then, to establish the variation in body temperature 
and water loss rates (EWLa and EWLi) according to popu- 
lation, sex and time interval, we used a (co)variance analy- 
sis with repeated measurement (AN(C)OVA-rm). Subse- 
quently, in both types of analysis (two-way and one-way 
AN(C)OVA and AN(C)OVA-rm), Log BM and Log SVL 
were incorporated as covariates to account for the effects of 
size, shape and body condition. Tukey HSD post hoc tests 


were performed between population pairs and sexes to ex- 
plore any significant differences. In the end, a partial corre- 
lation was performed between Tpref and EWLt at individ- 
ual level, taking SVL, BM, population and sex into account. 

All our analyses were computed in IBM SPSS Statis- 
tics version 22.0 (IBM Corporation 2011). Significance was 
evaluated at an a-value of 0.05. 


Results 


Our statistical analyses of snout-vent length (Table 1) re- 
vealed variation between populations (the Ait Bouguemez 
population being larger) and sexual dimorphism (females 
being larger), but no sexual dimorphism across popula- 
tions. The results for body mass were similar, except that 
there was no sexual dimorphism (two-way ANOVA test: 
for SVL, Populations a = 13.415, p < 0.001, Sex ae = 
6.941, p = 0.011, Populations x Sex LC 0.447, p = 0.642; 
for BM, Populations F, .,= 15-140, p < 0.001, Sex F, _,= 0.164, 
p = 0.687, Populations x Sex r= 1.177, p = 0.316). Simi- 
larly, our analysis of the body mass when SVL was used 
as a covariate showed a significant difference in robustness 
between populations and sexes, but no interaction between 
these factors (two-way ANCOVA test for BM, Populations 
a= 4.166, p = 0.021, Sex lg 9.523, p = 0.003, Popula- 
tions x Sex F, „= 1-035, p = 0.362). 

The analyses of two-way (co)variance of preferred body 
temperature and total water loss produced different re- 
sults for both parameters. The analyses of Tpref (Table 1) 
indicated no differences between the three populations, 
but did so between sexes (females having a lower Tpref) 
and no interaction Populations x Sex (two-way ANOVA; 
Populations ar 3.069, p = 0.055, Sex Bys 6.916, p = 
0.011, Populations x Sex F7 0.755; p = 0.475). Howev- 
er, after accounting for lizards’ size, the opposite pattern 
was found, with differences between populations, but nei- 
ther between sexes nor the interaction Populations x Sex 
(two-way ANCOVA with log BM and log SVL as covari- 
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ates, logBM F, _ = 0.001, p = 0.973, logSVL F, _ = 2.709, p = 
0.106, Populations F s2 = 5-975» P = 0.005, Sex Fo = 2.343 
p = 0.132, Populations x Sex F, ga = 0-415, p = 0. 662). Finally, 
ANOVA-rm (Fig. 2) indicated temporal variations in Tpref 
with different patterns for populations and sexes, which 
mostly disappeared after accounting for body size/shape 
(Supplementary Table S2). 

For EWLt (Table 1), the analyses of raw data revealed a 
difference between populations, but not between sexes or 
interaction Populations x Sex (two-way ANOVA; Popula- 
tions F, = 4.446, p = 0.016, Sex F, = 0.738, p = 0.394, 
Populations x Sex F, „7 2199, p = 0. 121). After accounting 
for size, no differences either between populations, sexes, 
or the interactions between them were found (two-way 
ANCOVA with log BM and log SVL as covariates, logBM 
F, s27 4255» P = 0.044, logSVL F, _, = 0.022, p = 0.882, Popu- 
lations F, ., = 0.400, p = 0.672, Sex F,_, = 0.307, p = 0.582, 
Populations x Sex F, ga = 1543, P = 0. 233). 

Through time, we found significant differences in all 
analyses with two-way ANOVA-rm (Time-Tpref, Time- 
EWLa and Time-EWLi, p<o.oo01). The analysis of Tpref 
produced significant differences in time between popula- 
tions as well as between sexes (ANOVA-rm; time x Pop- 
ulations F, osgo Z 1971, P = 0.007; time x sex S = 3.137, 
p=o. 001). In EWLa rates, we found differences in time 
between populations (ANOVA-rm; _ time*Populations 

325947 3:576; P = 0.019), but no differences between sexes 
(ANOVA-rm; time x sex F n97 0-807 p = 0. 413). In EWLi 
rates, no differences were detected in time either between 
populations or between sexes (ANOVA-rm; time x popu- 
lation F go 9 783, p = 0.606; time x sex F aig Os 617, p= 
0.634). Tie wike: no significant difference was found in the 
interaction time x populations x sex (p > 0.05) in Tpref, 
EWLi or EWLa (Figs 2-4, respectively). 

Finally, partial correlation analyses revealed a significant 
negative relationship between EWLt and logBM-logSVL 


Ait Bouguemez 


Preferred temperature (°C) 


Oukaimeden 


(n= 60; df =58,r en logBM-EWLt = -0.518, p < 0.001; df = 
58, T paria OSV L- -EWLt = -o. 369, p = 0.004), with all pop- 
ulations following the same trend. In contrast, no signifi- 
cant association of Tpref with either logBM or logSVL was 
observed (n = 60; df = 58, T ma logBM-Tpref = -0.017, p = 
0.895; df = 58, r ed logSVL-Tpref = -0.182, p = 0.165). No 
significant association was observed even after accounting 
for logBM, logSVL, populations and sex (n = 60; df = 54, 
LP sig) EWLt-Tpref = 0.146, p = 0.282). 


Discussion 


Our ecophysiological comparisons across three widely di- 
vergent lineages within the Atlantolacerta andreanskyi spe- 
cies complex, encompassing more than 7 Mya of independ- 
ent evolution, is in line with previous findings suggesting 
that thermal physiology is more rigid than hydric physiol- 
ogy is along the phylogenetic axis in Lacertidae (GaRcIA- 
Porta et al. 2019). However, they also provide evidence of 
phenotypic plasticity, since thermoregulation in gravid fe- 
males was biased towards lower temperatures (CARRETERO 
et al. 2005) and mainly because variations in body size also 
had repercussions on the lizards’ hydric physiology (FER- 
REIRA et al. 2016). 

As in other lacertids from cold climates, sexual dimor- 
phism in SVL of A. andreanskyi indicates that females 
are larger, likely due to stronger fecundity selection rath- 
er than sexual selection under these conditions (BRANA 
1996, ROITBERG et al. 2015). Moreover, in populations with 
shorter annual activity periods, growth takes precedence 
over reproduction and lizards tend to delay sexual matu- 
rity in favour of attaining larger sizes first (HORVÁTHOVÁ 
et al. 2013, ROITBERG et al. 2020). Whereas robustness of 
male lacertids compared to females is considered a result of 
sexual selection (KALIONTZOPOULOU et al. 2008), robust- 
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Time (Hours) 


Figure 2. Daily variations in the preferred body temperatures (Tpref) of Atlantolacerta andreanskyi (males left, females right) between 
the three studied populations and sexes. Displayed are mean values and standard deviations. 
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ness of both sexes rising with altitude is probably the result 
of allometric growth (KALIONTZOPOULOU et al. 2010a, b). 
Despite their considerable genetic separation, preferred 
body temperatures were similar across the studied popula- 
tions of A. andreanskyi. This conservativeness is not sur- 
prising and has been observed along other wider phylo- 
genetic distances in this lizard family (GARCIA-PORTA et 
al. 2019). However, the preference of lower temperatures 
by (gravid) females indicates there is phenotypic plasticity 
also for this trait. In Mountain Lizards with short activ- 
ity seasons, female phenology tends to be accentuated and 
synchronized (CARRETERO et al. 2005, 2006), suggesting 
the gestation status to have been similar across examined 


Ait Bouguemez 


Oukaimeden 


Instantaneous evaporative water loss rate (%) 


andreanskyi to climate change 


populations and individuals, and mirrored by similar ther- 
mal shifts across populations. While an effect of thermal 
stress on advanced embryos cannot be ruled out (MADER- 
SON & BELLAIRS 1962), the fact that the statistical removal 
of size/shape effects made the shift disappear suggests oth- 
er causative factors, namely, hydric physiology. 

Like all lacertids, all populations of A. andreanskyi 
tended to raise their Tpref in the early morning when en- 
vironmental temperatures in the wild would be low, that 
is, when the thermal environment is more unfavourable 
and thermoregulation becomes a priority relative to oth- 
er physiological requirements (RISMILLER & HELDMAIER 
1982, TOSINI & AVERY 1994, CARRETERO et al. 2005). In- 
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Figure 3. Daily variations in the instantaneous water loss rates (EWLi) of Atlantolacerta andreanskyi (males left, females right) between 
the three studied populations and sexes. Displayed are mean values and standard deviations. 
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Figure 4. Daily variations in the accumulated water loss rates (EWLa) of Atlantolacerta andreanskyi (males left, females right) between 
the three studied populations and sexes. Displayed are mean values and standard deviations. 
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terestingly, this was more accentuated in males than in fe- 
males. However, the fact that Tpref values in the late after- 
noon tended to decline rather than increase suggests that 
water loss was constraining thermoregulation towards the 
end of the day (S’KH1Fa et al. 2020a). Once more, these 
shifts were dependent on lizard size and shape and mostly 
disappeared after accounting for SVL and mass. In a study 
comparing lacertid species representing different sizes, cli- 
matic regions and clades, thermal and water ecophysiology 
showed signs of a trade-off, but patterns varied between 
species and were independent of size/form (FERREIRA et 
al. 2016). Here, the trade-off was also evident (lizards with 
high Tpref had lower EWL and vice versa), but all three 
populations followed a common trend and the trade-off 
was clearly size/form-dependent, pointing again to con- 
servatism of fundamental niche occupancy. 

A strong population signal was recovered for EWL in 
A. andreanskyi, namely, a decrease with increasing alti- 
tude. These results were expectable since SANNOLO et al. 
(2020) already demonstrated flexibility in hydroregula- 
tion by the lacertid Psammodromus algirus from different 
altitude across a small geographic gradient. While the ef- 
fect of local acclimatization (e.g., to different wind inten- 
sities) and the separate evolutionary histories of different 
lineages cannot be ignored (Dm’EL et al. 1997, Cox & Cox 
2015), the simplest explanation in this case is that these dif- 
ferences in EWL result from body size increasing with in- 
creasing altitude, which decreased the surface/volume ra- 
tio (SCHMIDT-NIELSEN 1984). Accordingly, the largest liz- 
ards from Ait Bouguemez exhibited considerably more 
resistance to water loss than the smallest ones from Jbel 
Aoulime, while those from Oukaimeden had intermediate 
values. This is confirmed by the lack of differences in EWL 
when SVL and body mass were accounted for. Remarkably, 
sex had no effect per se on the results, although the larg- 
er sizes and more voluminous abdomens of females may 
have minimised their water loss as a function of gravid- 
ity (BONNET et al. 2017). This suggests that thermal ecolo- 
gy is taking precedence over hydric ecology regarding the 
embryo development inside the mother. One possible ex- 
planation would be the reduced oxygen pressure at very 
high altitudes (GANGLoFF et al. 2019). Certainly, after ovi- 
position and then during the juvenile stage it is clear that 
hydroregulation will increase in importance. Evidence for 
this is the finding of communal nests and juvenile aggrega- 
tions, which suggest that suitable nesting sites are limited 
(Marco et al. 2004, S KHIFA et al. 2020b). 

Behavioural thermoregulation is understood to be a key 
mechanism for buffering the physiological effects of cli- 
mate change, but the classic overreliance on thermal test- 
ing without accounting for water availability have limited 
the development of generalised models (ROZEN-RECHELS 
et al. 2020). Our data, combined with other recent evalua- 
tions of EWL in lizards (e.g., BELASEN et al. 2017), further 
support the general pattern that hydric ecology is more 
flexible than thermal ecology, and therefore is the key com- 
ponent to understanding the impact of climate change on 
ectothermic vertebrates. The development of worldwide 
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databases for EWL rates in squamates (e.g., LE GALLIARD 
et al. 2021) may help further address this issue, although 
undersampling is evident in various geographic regions, 
including Africa. Rather than focussing primarily on high- 
er temperatures and available thermal niches, a more nu- 
anced approach including water availability and microhab- 
itat will be needed to better model the impact of climate 
change on lizards in general. 

Our findings also have important repercussions on con- 
servation measures for our target species. Atlantolacerta 
andreanskyi is the only representative of its genus (ARNOLD 
et al. 2007), sister taxon to all other Eremiadini in the lac- 
ertid phylogeny (GarciA-Porta et al. 2019), and consid- 
ered as Near Threatened in the Red List of IUCN due to its 
reduced and fragmented range and high habitat specific- 
ity (GENIEZ 2006). The ecophysiological evidence reported 
here is also consistent with correlative ecological models, 
underscoring the importance of precipitation for this spe- 
cies’ occurrence and predicting a considerable range retrac- 
tion under the climate change scenarios (MARTiNEZ-FEIRIA 
et al., 2013). However, this considers the complex as a whole 
- if the genetic lineages were formally described as distinct 
species, some would likely need to be reclassified as En- 
dangered due to their small known ranges, especially giv- 
en that the range of the whole complex is already less than 
20,000 km’. Because it lives in harsh environments at high 
altitudes in the Atlas Mountains, its vulnerability to global 
warming is expected to be very high (SINERVo et al. 2010). 
Niche conservatism in thermal ecophysiology within this 
species complex reduces the extent of expectations for fast 
evolutionary responses to fast environmental shifts. Plastic- 
ity in hydric ecophysiology between populations and sexes 
might provide some optimism, but this was mainly mediat- 
ed by body size, and the smallest lizards of the populations 
at lower altitudes would in fact be the ones that are most 
vulnerable to aridification (HIJMANS et al. 2005, S KHIFA 
et al. 2020a). Furthermore, thermal constraints imposed 
on embryos duringgestation, as well as the hydric needs of 
eggs and hatchlings in an environment with limited hydric 
resources may compromise reproduction and pose an addi- 
tional threat. While little can be done locally to reverse the 
global trends in climate change, conserving natural habi- 
tats, and in particular the vegetation coverage and incor- 
porated water bodies, seems the best management strategy 
at a local level to prevent fast range retractions, population 
declines and, ultimately, extinction (GENIEZ 2006). Part of 
the range of this Moroccan endemite lies within a protected 
area, the Toubkal National Park. However, this would only 
protect one (Oukaimeden) of the several identified genetic 
lineages. For the remaining populations, where no conser- 
vation measures are currently in place, habitat amelioration 
focusing on vegetation and habitats that offer moisture is 
strongly recommended. The IUCN Red List does not cur- 
rently report any specific threats, although it notes that the 
species occurs in areas accessible to hikers. Other than in 
the national park, we suggest that very few hikers traverse 
the places where these populations occur. Rather, given the 
importance of source of moisture in habitats, efforts could 
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be focussed on reducing overgrazing in these fragile lo- 
calities. Our findings are expected to improve our under- 
standing of the potential repercussions of climate change 
on the conservation status of this Moroccan endemic (GE- 
NIEZ 2006; MARTÍNEZ-FEIRÍA et al., 2013), supply scientific 
evidence for its management, and provide further data to 
develop generalized models for assessing the impact of cli- 
mate change on ectothermic species. 
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